Approaching the ground state of the kagome antiferromagnet 
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Yo.5Cao.5BaCc>407 contains kagome layers of Co ions, whose spins are strongly coupled according 
to a Curie-Weiss temperature of -2200 K. At low temperatures, T = 1.2 K, our diffuse neutron 
scattering study with polarization analysis reveals characteristic spin correlations close to a predicted 
two-dimensional coplanar ground state with staggered chirality. The absence of three dimensional 
long-range AF order proves negligible coupling between the kagome layers. The scattering intensities 
are consistent with high spin S = 3/2 states of Co 2+ in the kagome layers and low spin S = states 
for Co 3+ ions at interlayer sites. Our observations agree with previous Monte Carlo simulations 
indicating a ground state of only short range chiral order. 

PACS numbers: 75.25.+Z, 75.40.Cx, 75.50.Ee 
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In low dimensions, antiferromagnetic (AF) order is 
suppressed at finite temperatures [1| and geometrical 
frustration raises the complexity of ground states with 
finite entropy and non-collinear, chiral spin structures 
3] . Spin chirality is related to unique universality classes 
sjlj] and has relevance for the spin excitations in copper- 
oxide superconductors d, @. A challenge for theoreti- 
cal understanding for decades is the utmost frustration 
among AF coupled spins on the two-dimensional kagome 
lattice 0- 

Considering the configurational energy of the classical 
Heisenberg AF, 
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for AF interaction to only nearest neighbors, J(|n.|) > 0, 
the ground state is highly degenerate. There are two 
competing ordered coplanar spin structures (see Fig. 1), 
with relative spin orientations of 120 degrees, which are 
either of uniform or staggered chirality. We use the con- 
vention that the chirality is positive when the spins ro- 
tate in steps of 120 degrees as one goes around a triangle. 
Apart from the degeneracy of a common in-plane spin ro- 
tation, there are further degeneracies related to possible 
local disorder, as shown in Fig. 1, Within its configura- 
tional space, thermal fluctuations 10, 11] as well as quan- 
tum fluctuations || select the V3 x -\/3 structure. The 
entropical selection resolving a ground state degeneracy 
is a mechanism that has been named order by disorder 
[l2|. The classical ground state has still a macroscopic 
manifold, which could be removed by quantum fluctu- 
ations 0. On the other hand, there are predictions of 
a disordered ground state, a quantum spin liquid with 
singlet formation for the S = 1/2 kagome AF |l3| . 

The topology of the classical ground state degeneracies 
has been analysed in more detail in Monte Carlo (MC) 
simulations of the 2D Heisenberg AF by Reimers and 
Berlinsky [JjJ . An intriguing result of these simulations 
is (an indication of) a divergence of the spin-correlation 




FIG. 1: (color online). Infinitely degenerate AF ground states 
on the kagome lattice, (a) the go-structure with uniform chi- 
rality; (b) the staggered chiral structure, which has a larger 
unit cell by y/3 x Chirality is denoted by signs inside the 
triangles. The dashed ellipses represent thermal fluctuations 
of zero energy: (a) common spin rotations in an infinite chain; 
(b) the weathervane defect, a common local spin rotation on 
the hexagon around the axis of the spins on the edges of the 
unit cell. 



function for T — > K towards an ordered ground state 
resembling the v3 x y3 structure, whereas chiral long- 
range order appears to be prohibited by disorder due to 
the local common zero energy spin rotations as depicted 
in Fig. 1. 

Experimental studies of the 2D kag ome AF are difficult 
and rare, although pyrochlores |14j . magnetoplumbitcs 
[l5j . and jarosites contain layers of the kagome 

structure. For weak interlayer coupling, these candidates 
combine all the desired ingredients: low dimensional- 
ity, low coordination, and strong frustration. A neutron 
scattering study on magnetoplumbite SrCrs-zGa^^Oig 
[lfij provided the first experimental evidence of short- 
range spin correlations down to 1.5 K and indicated 
a preference for the V3 x structure. The jarosite 
(D30)Fe3(S04)2(OD)6 stays disordered even at low T ac- 
cording to diffuse neutron scattering data (obtained with 
polarization analysis), but the ordering wave vector could 
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FIG. 2: (color online). The average crystal structure of 
Yo.5Cao.5BaCo407. Thicker lines between the Co2 atoms 
highlight the kagome layers. All Co atoms are tetrahedrally 
coordinated by oxygen, as shown in the polyhedral represen- 
tation (left), emphasizing the three-fold symmetry axes, going 
through Col (blue), which allow trimerization of Co2 in the 
kagome substructure. 



not be clearly identified. Another close realization of the 
kagome AF is the jarosite KFe3(S04)2(OD) 6 , although 
the interference of 3D Neel order below 65 K complicates 
comparisons to the ground state properties of the kagome 
AF. Recent single crystal studies show spin correlations 
towards the go-structure |l7| . contrary to theoretical ex- 
pectations for the kagome AF with only nearest neighbor 
interactions. Indeed, small AF interactions to next near- 
est neighbors, which have been deduced from the low en- 
ergy spin wave excitations [lij , stabilize the qo-structure 
instead of an order by disorder mechanism. More inter- 
estingly, the low energy spin wave excitations could be 
related to the "zero energy modes" as shown in Fig. 1 (a) 



Here, we present results of diffuse neutron scatter- 
ing with polarization analysis on a new compound 
Yo.5Cao.5BaCo407, a realization of the kagome AF with 
only nearest neighbor interactions, which allows an un- 
precedented approach to its ground state properties. 

The average crystal structure of Y .5Cao.5BaCo407 is 
known [l^. E(i| (Fig. 2) and a mixed valence of Co 2+ 
and Co 3+ is a consequence of stoichiometry. Co occu- 
pies two different tetrahedral sites in the structure: 75% 
of all Co (Co2) constitute the kagome layers and the re- 
maining 25% fills the Col sites, which are on threefold 
axes between the kagome layers. Within this symmetry, 
a trimerization in the kagome layers is possible and dif- 
ferent Co2-Co2 bond lengths are in fact observed at 300 
K [23. Because of shorter Co-0 distances, the valence 
of Col must be higher than that of Co2, indicating a 
preference for Co 3+ (3d 6 configuration). The bonding 
situation and spin state of Co 3+ seem to be responsible 
for two dimensional magnetic behavior: (i) geometrically, 
the tetrahedral coordination causes a high frustration for 



any possible AF exchange between the kagome layers; (ii) 
the possibility of a low spin state of Co 3+ at the inter- 
layer sites, as Co 3+ is situated closer to a trigonal plane 
changing the point-group symmetry from to C% v . Ac- 
cording to Hund's rules, an ideal tetrahedral coordination 
with weak crystal field splitting (oxygen ligands) results 
in high spin states, on the other hand low spin states 
have been observed for C^ v symmetry [^, [2^ • 

Polycrystalline Yo.5Cao.sBaCo407 samples for the dif- 
fuse neutron scattering experiments have been synthe- 
sized by a solid-state reaction in air [2(J. Under the 
applied synthesis conditions, the oxygen stoichiometry 
should be ideal as verified by iodometric titration on 
the parent compound YBaCo4C>7 j^. The scatter- 
ing experiments were performed on the DNS instrument 
in Julich. In contrast to most previous investigations 
13 US- Ei j the magnetic scattering was separated from 
nuclear scattering by polarization analysis |16t |25j. A 
fortiori, the separated high background due to nuclear 
spin-incoherent scattering from Co serves as an intrinsic 
calibration of the paramagnetic cross-section; a proce- 
dure that avoids usual systematic errors and corrections. 
A strong diffuse peak evolves upon cooling and its asym- 
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FIG. 3: (color online). Separation of scattering contributions 
and temperature dependence of the separated diffuse mag- 
netic scattering; data are calibrated by the separated spin- 
incoherent scattering of Co (0.382 b/sr); the separated nu- 
clear coherent scattering, I nc has been rescaled and shifted 
(bottom part). 



metric shape indicates low dimensional spin correlations 
(Fig. 3). Remarkably, no long-range magnetic order ap- 
pears even at very low T. The magnetic scattering can be 
described in terms of the Fourier transform of the spin 
pair-correlations, appl ying orientational averaging for the 
polycrystalline sample |2J] . For a direct comparison with 
the MC results, we use a simplified scattering expression 
|ll| . which neglects a possible interrelationship between 



3 



spin correlation and spin direction: 
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Here, the product of the classical electron radius 
(e 2 /m e c 2 ) = 0.282 x 10~ 12 cm and the neutron gyromag- 
netic ratio = -1.91 defines the magnetic scattering length 
i'q = —0.54 x 10~ 12 cm; m e is the electron mass, e is 
its charge, and c is the velocity of light. S represents 
the spin quantum number of the scattering ion, Fq de- 
notes the magnetic scattering form- factor of the single ion 
Co 2+ . For the ideal paramagnet, T> J, the scattering 
is purely clastic and proportional to the self-correlation 
S(S + 1). In general, one has to consider the expectation 
value of the time-dependent self-correlation function of 
a spin precessing in the strong local field of its neigh- 
bors. For T< J and integrating over only small energy 
transfers in diffraction, the expectation value reduces ef- 
fectively to that of the scalar product of the ordered mo- 
ments {Sf). 

The magnetic scattering intensity allows us to esti- 
mate the average Co spin moment. A fit to the 1.2 K 
data yields a forward scattering of 0.25(2) b/sr-Co for 
the self-correlation term. The best quantitative agree- 
ment is found if we assume that Co 2+ ions (5/8 of all 
Co) are in high spin state S = 3/2 and that Co 3+ ions 
(3/8 of all Co) are in low spin state 5 = yielding 0.273 
b/sr-Co-ion. It is clear that S = at the Col sites ratio- 
nalizes the absence of any significant interlayer coupling 
and, hence, the ultimate suppression of 3D long-range 
order. 

The diffuse magnetic scattering probes the wave- vector 
dependent susceptibility x{Q) = T^ t {Sq- S r ) slr ^ r , where 
S are spins of unit length. Within the approximation 
made for Eq. (2), the 1.2 K data are compared with the 
MC results in Fig. 4. Correlations between the direc- 
tions of the spins and their distance vector would require 
higher order corrections neglected in Eq. (2). In case of 
the \/3 x \f?> structure, as shown in Fig. 1(b), we ana- 
lyzed that this would result in significant intensities for 
< Qr < 1, for nearest neighbor distances |ri| = a/2 
and a unit cell parameter a = 6.30 A, which is in con- 
tradiction to our observations. According to Eq. (1), the 
Hamiltonian is degenerate with respect to a common spin 
rotation in the kagome plane, and we conclude that this 
degeneracy is essentially preserved. 

Furthermore, we can rule out spin correlations to- 
wards the qo-structure, whose ordering wave vectors are 
a*(h,k); the strongest peak would be at Q = |qo| = 

|a*(l, 0)| = 1.152 A" 1 , where a* = (2 7 r/a)(2/ v / 3). In- 
stead of this, the observed peak coincides with Q — 
|2q^| = |2a*(l/3, 1/3) | = 1.330 A~\ the ordering wave 
vector of the V3 x y3 structure. 

The data has been analyzed in terms of spin correla- 
tions by a Fourier analysis using a linear least squares 



refinement. The results are shown in the inset of Fig.0] 
where the spin correlations have been normalized to the 
correlation function of the ideal x y/3 structure: 



G v^W = (So ' SVVcosq^g 



(3) 



which equals to unity when the spins are in the v3 x \/3 
configuration. A limited number of parameters could be 
determined, which clearly show the preference for build- 
ing up the -\/3 x -\/3 structure. 
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FIG. 4: (color online). Susceptibility x(Q) from 1.2 K data 
versus modulus of wave vector Q; the thick (cyan) line rep- 
resents MC data from Ref. [ll|. for a classical Heisenberg 
AF at T/J = 0.002; dashed line represents a fit of the spin 
correlation by Fourier analysis, whose results, as normalized 



to correlations in the ideal structure, 
inset. 
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An interesting question is whether our data can pro- 
vide any indication for the proposed weathervane defects 
(cf. Fig. 1(b)) that seem to disturb the chiral long range 
order of the ground state in the previous MC simula- 
tions. Therefore, one may compare a model of the ideal 
\/3 x \/3 structure to a defect model, and convolute the 
correlations with an adjustable asymptotic decay func- 
tion. A simple defect model is to place the weathervane 
defects on two of the three sublattices. Note that the 
configurational average of this model is equivalent to a 
collinear structure model, and that the correlations are 
reduced by a factor 2 with respect the ideal structure ( 
(So ■ S r ) = 1 or —1/2) for distances exceeding the spac- 
ings that occur within the unit cell. An algebraic decay 
of the correlations has been found in MC simulations [ll| , 
which relates the spin correlations to a percolation phe- 
nomenon. However, the present data are affected by in- 
strumental resolution (see Fig. 3 separated nuclear scat- 
tering) such that we cannot reliably distinguish between 
an exponential and algebraic decay. Models based either 
on the ideal structure or on the weathervane defect are 
in close agreement with the data. The exception is that 
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both models produce a small additional peak at q^, an 
integral intensity that is one order of magnitude weaker 
than the peak at 2q^. The discrepancy is less signifi- 
cant, approximately half the size, if the weathervane de- 
fects are included. Apparently, the spin structure is even 
less ordered and a more appropriate model would include 
that weathervane defects are distributed on all sublat- 
tices. Therefore, MC simulations are adequate to sample 
the configurational space. Indeed, there is a qualitatively 
astonishing agreement between our data and the MC re- 
sults for the classical Heisenberg AF. In particular, 
MC simulations reproduces the absence of the peak at 
O/g. The disorder has been related to a structure of chiral 
domains with anti-phase boundaries. According to the 
MC data, the correlation length of chiral order parame- 
ter stays short ranged for T — ► K , while the spin corre- 
lation length and the peak at 2g^g diverge in this limit. 
It is noteworthy that this observed peak, as compared 
to the MC simulation, is significantly more pronounced. 
Hence, the ground state is approached much closer in 
our experiment, although temperature and Curie- Weiss 
constant of approximately — 2200 K w 4 J |2(| matches 
closely the conditions of the simulation, T/J= 0.002. In 
fact, the peak is affected by the instrumental resolution 
and the true peak height is already substantially reduced. 

The weathervane defects (see Fig. 1(b)) are local collec- 
tive spin excitations of zero energy according to Eq. (1) 
that preserve the sum of spins and give rise to residual 
entropy in analogy to ice [2tJ and spin- ice systems [l4j| . 
Local collective spin excitations have also been observed 
in inelastic neutron scattering from cubic spinel ZnCr2 04 
|26| . Very recently, "chain- like" excitations (from the go- 
structure, see Fig. 1(a)) are lifted to finite energies by 
Dzyaloshinskii-Moriya (DM) interaction as been reported 
from inelastic neutron scattering [J^. We performed 
additional time-of-flight experiments verifying that the 
scattering at 1.2 K is essentially elastic, and proving that 
the spin correlations are frozen within the energy reso- 
lution of 0.5 meV. At higher temperatures, the inelas- 
tic scattering appears gapless, indicating the absence of 
anisotropics related to DM interaction, and exhibits a 
quasi-elastic broadening similar to the observations on 
SrCrg^Ga^Org O. 

The strong AF exchange in the plane causes a dynamic 
response at high energies, which suggests that the low- 
energetic relaxations correspond to thermal spin fluctua- 
tions out of the kagome plane. Theoretically, the complex 
free energy landscape of the kagome AF will induce very 
slow spin dynamics similar to ordinary spin glasses |2Sj 
and may explain observed field cooling effects at 370 K 

To conclude, the present neutron study with polariza- 
tion analysis on a new compound, realizing the 5 = 3/2 
kagome AF, agrees in surprising detail with the complex- 
ity of the classical ground state for only nearest neighbor 
interactions as has been found in previous MC simula- 



tions. The remarkable survival of only short-range 2D 
spin correlations at low temperatures is due to severe 
geometrical frustration by tetrahedral coordination com- 
bined with high/low spin states of Co ions in alternating 
layers. Site disorder with spin dilution, albeit beyond 
the site percolation threshold [2^| , results from assuming 
ionic bonding and from exchange of moments through 
electron hopping. Covalency effects might restore the 
predicted behavior in the experiment, which implies a 
long-range ordered ground state with only short-range 
chiral order. 

We thank J.N. Reimers and A.J. Berlinsky for provid- 
ing us with their original MC data. 
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